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C^ . PACS. 71.20. Be - Transition metals and alloys. 

PACS. 71.27.+a- Strongly correlated electron systems; heavy fermions. 
PACS. 75.20.Hr - Local moments in compounds and alloys. 

o 

Q . 

Abstract. - We present results of all-electron electronic structure calculations for the recently 
discovered d electron heavy fermion compound LiV2 04. The augmented spherical wave cal- 
culations are based on density functional theory within the local density approximation. The 
^ ' electronic properties near the Fermi energy originate almost exclusively from V 3d t2g states, 

\l I which fall into two equally occupied subbands: While cr-type metal-metal bonding leads to 

0^ , rather broad bands, small 7r-type p-d overlap causes a narrow peak at Ep. Without the 

geometric frustration inherent in the crystal structure, spin-polarized calculations reveal an 
>0 ' antiferromagnetic ground state and ferromagnetic order at slightly higher energy. Since direct 

^^ I d-d exchange interaction plays only a minor role, ordering of the localized vanadium moments 

0^ , can be attributed exclusively to a rather weak superexchange interaction. With the magnetic 

0^ ■ order suppressed by the geometric frustration, the remaining spin fluctuations suggest an ex- 

planation of the low temperature behaviour of the speciflc heat. 

a 

O ■ The magnetic porperties of the metallic spinel structure transition metal oxide LiV2 04 have 

remained an unresolved issue for some time [nl 0, pi B . Apart from a temperature independent 

contribution, the susceptibility was found to follow a Curie- Weiss law due to local V magnetic 

k><( \ moments O, El . The negative Weiss temperature of -63 K indicates antiferromagnetic exchange 

5_j ' interactions between the vanadium spins H. However, a transition to long-range magnetic 

C^ \ order above 1.8 K could be excluded from both susceptibility and ^Li NMR measurements 

[0, o| . This apparent contradiction was attributed to possible geometric frustration inherent in 

the vanadium sublattice of the normal spinel structure with only nearest neighbour interaction, 

which might prohibit long-range antiferromagnetic order H, 0| . 

Recently, Kondo et al. reported a crossover from localized moment to heavy Fermi liquid 
behaviour with a Kondo temperature of Tk ~ 28 K and an electronic specific heat coefficient 
of 7 « 0.42J/molK^ at 1 K p, M. Thus for the first time heavy fermion (HF) behaviour, 
characteristic of / electron systems, has been observed in a d electron material. Below 30 
K, the resistivity exhibited a pronounced smooth downturn consistent with Kondo lattice 
behaviour [H. The nearly temperature independent susceptibility and the ^Li Knight shift 
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in the low temperature range [§| ^, Ig] were interpreted as due to the disappearance of the V 
local moments. Neither static magnetic order nor superconductivity was observed above 0.02 
K g. Yet, low temperature /iSR data displayed anomalies indicative of static spin freezing 
below 0.8 K Ig, ^. The specific heat shows an upturn between 0.8 and 0.5 K, which might 
point to a T'^ In T dependence characteristic of large spin fluctuations [Q . Low temperature 
^Li and ^^V NMR data likewise showed deviations from uniform heavy Fermi liquid behaviour 
[ [ll[ [l^ , |l3|. From neutron diffraction experiments, structural transitions between 4 and 295 K 
could be excluded |l4]. In addition, these measurements revealed a large electronic Griineisen 
parameter comparable to that of / electron HF compounds, which indicates strong coupling 
of the electronic and lattice degrees of freedom H, M. In contrast to the / electron systems, 
however, where HF behaviour as well as long-range antiferromagnetic order result from the 
hybridization of localized and bandlike orbitals, the origin of the electron-lattice coupling and 
the mechanisms leading to HF formation in the d band metal LiV2 04 are as yet unclear [Q. 

In this Letter, we present for the first time results of all-electron first principles calculations 
for LiV2 04. Our extensive study is based on density functional theory (DFT) within the local 
density approximation (LDA) ||l^. As a calculational scheme we used the scalar-relativistic 
augmented spherical wave (ASW) method ||l6[. As a result, we find (i) a superposition of 
narrow and broad bands at the Fermi energy. Allowing for magnetic ordering of the vanadium 
moments, we obtain (ii) an antiferromagnetic ground state resulting from weak superexchange 
coupling, which is accompanied by ferromagnetic order at a slightly elevated energy. Due to 
the geometric frustration only large spin fluctuations remain, while the reduced overlap of the 
magnetic orbitals points to a possible localization transition. 

LiV2 04 crystallizes in the normal spinel structure, which is based on an fee lattice with 
space group FdZm (O^) H, |l^. We used the values of a = 8.22672 A for the lattice constant 
and XQ — 0.26111 for the oxygen parameter as resulting from the 4 K neutron diffraction 
data llj]. In the normal spinel structure, the oxygen atoms form an almost perfect cubic 
close-packed sublattice with tetrahedral and octahedral interstices. The lithium atoms occupy 
1/8 of the tetrahedral sites whereas the vanadium atoms occupy half of the octahedral sites 
and form two interpenetrating fee sublattices of corner sharing tetrahedra with a V-V bond 
length of 2.9 A B. Due to cation ordering the crystal field at the transition metal sites is 
trigonal rather than cubic. The deviation of the oxygen xq parameter from the ideal value 0.25 
causes a shift of 0.158 A in (111) direction, which enlarges/compresses the available volume at 
the tetrahedral/octahedral sites; the V-0 distance is decreased by 4.6%. At the same time, 
the V-O-V bond angle changes to 95.2°. 

Fig. 1^ displays the dominant partial densities of states (DOS). States not included play only 
a negligible role in the energy interval shown. Three groups of bands are identified: While 
bands in the energy range from -8 to -3.3 eV originate mainly from O 2p states and have only 
a small admixture from V id states, the upper two groups of bands, which extend from -1.3 to 
0.8 eV and 2.4 to 3.4 eV, are predominantly derived from the V "id states. Although additional 
O 2p contributions are apparent in this energy range, p-d hybridization is much reduced as 
compared to other early transition metal oxides |18]. As concerns the occupied states our 
findings are in good agreement with the photoemission data |l9]. Octahedral crystal field 
splitting, as expected from the nearly cubic coordination of the V atoms by the oxygen atoms, 
leads to clear energetical separation of the Zd i2g and Cg groups of bands as is visible in 
fig. H. Whereas the former states appear exclusively around Ep, the Cg states dominate at 
higher energies. Contributions of the V id states to the oxygen derived bands originate almost 
exclusively from the eg states, which, forming <j bonds, experience a larger overlap with the 
O 2p states. In contrast, the t2g orbitals, which give rise to tt bonds, yield only a negligible 
contribution in this energy range. In addition to these rather weak p-d bonds the t2g states 
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Fig. 1. - Partial densities of states (DOS) of LiV204 per unit cell. 



experience, however, strong cr-type overlap with the t2g orbitals at neighbouring vanadium 
sites of the fee sublattice. Hence, these d states take part in two different types of bonding, 
namely cr-type V~V- and 7r-type V-0-bonding, which leads to two different band widths. 
Nevertheless, since both the metal-metal and the metal-oxygen bonding arc mediated by the 
same orbitals, a simple analysis of the partial DOS would not allow to distinguish the different 
roles played by the t2g orbitals. 

We solved this problem by complementing the calculations for the observed crystal structure 
(xQ = 0.26111) with another set, where the oxygen positional parameter was hypothetically 
adjusted to its ideal value of xq = 0.25. Since this leaves the vanadium-vanadium bonding 
unaffected, the changes of the electronic structure due to the oxygen shift allow to make the 
aforementioned distinction. The results arc shown in fig. 2, where we display the partial DOS 
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Fig. 2. - Partial densities of states of LiV204 for the (a) observed (xq = 0.26111) and (b) ideal 
(xo ~ 0.25) structure. 
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for the V 3d dominated bands. We used the trigonal representation of the fcc lattice and 
plotted the V d^z^-r^, dxz + dyz and dx2_y2 + dxy partial DOS. While the d^^^-r^ orbitals, 
which correspond to the aig orbitals in a trigonal crystal field, point along the (111) direction 
and are of pure t2g character, the other four orbitals comprise a mixture of t2g and eg states. 
However, the latter states are dominated by the dxz and dyz orbitals, which are almost empty. 

For xo — 0.26111 we observe a clear separation of the t2g states into an Ri2.1 eV wide 
group formed by the dx2^y2 /dxy states and another group due to the much narrower d2,z^-r^ 
bands, which, apart from the peak at 0.7 eV, extend over an energy range of only 0.6 eV. 
On the other hand, for the idealized structure this difference of the band widths is removed 
due to the striking change of the d^^z^-r^ partial DOS: The peak at «0.35 eV is replaced by a 
nearly rectangular shaped DOS of wl.2 eV width. While the O 2p DOS closely follows these 
changes, the partial DOS due to the other d orbitals are almost unaffected by the oxygen 
atom shift. From inspection of the band structures the drastic changes of the d^z'^-r'^ states 
can be attributed almost exclusively to the striking downshift of a single band by up to 0.45 
eV. A similar situation has been found for LiTi204 ||2y, |2J. While in the real structure, this 
band Hes completely above Ep and causes the peak of the d^z^-r^ partial DOS at 0.35 eV, it 
is stabilized in the ideal structure by the widening of the triangular octahedral faces, which 
point towards the (111) direction. As a consequence, the electron count of these orbitals is 
slightly increased at the expense of the dx2-y2 and dxy electrons. 

In conclusion, from their different response to the oxygen atom shifts we are able to clearly 
identify two groups of t2g bands: The broad dx2-y2/dxy states originate from strong a-type 
V-V-overlap and are thus essentially unaffected by the change of the oxygen xq parameter. 
In contrast, the much narrower d^z^-r^ bands are involved in weak 7r-type V-0-overlap and 
depend sensitively on the oxygen position. 

Motivated by the observed Curie- Weiss behaviour and disregarding for the time being 
possible geometric frustration effects in the spinel structure, we performed additional spin- 
polarized calculations. Both ferro- and antiferromagnetic order with a modulation of spins 
parallel to one of the cubic axes (q = (0, 0, Air/a)) were allowed for and, again, we considered 
both the observed and the idealized crystal structure. In all cases we obtained magnetic order 
with vanadium magnetic moments and total energy lowerings as summarized in table |. The 



Table I. - Magnetic energy gain and vanadium magnetic moments. 





Xo 


= 0.26111 


Xo 


= 0.25 




ferro 


af 


ferro 


af 


Etot mRyd/f.u. 


-1.431 


-3.655 


-2.701 


-1.170 


mv [mb] 


0.590 


0.877 


1.108 


0.768 



latter are given relative to the respective non-spinpolarized solution. The vanadium magnetic 
moments result from almost equal contributions from the dx2_y2 /dxy and 1^322 _^2 orbitals; 
moments on all other atomic spheres are smaller than 0.025 /ib. The main results are: For the 
antiferromagnetic states the rather small change of the ordered moment with respect to shifting 
the oxygen xq parameter is contrasted with a considerable modification of the magnetic energy 
gain. In contrast, we witness a suppression of both the ferromagnetic moment and energy gain 
on going from the idealized to the observed structure. A detailed analysis revealed that this 
loss of magnetic moment affects all the t2g states to a similar degree. 

The aforementioned behaviour arises from two different mechanisms: (i) The isotropic 
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compression of the VOg octahcdra in the observed structure as indicated by the 4.6% decrease 
of the V-0 bond length, causes destabihzation of both the ferromagnetic moment and order as 
compared to the ideahzed structure. Furthermore, since the magnetic moments carried by the 
dx2 _y2 / dxy and d^^2_j.2 orbitals are affected to a similar degree by the oxygen atom shift, we 
conclude that the polarization of the former orbitals is due to Hund's rule coupling rather than 
direct exchange interaction across the cr-type V-V bonds, (ii) The deviation of the V-O-V 
bond angle from 90° for xq = 0.26111 supports antiferromagnetic superexchange interaction 
of the local vanadium magnetic moments via the occupied O 2p orbitals, while changing their 
size only little. The crossover from ferro- to antiferromagnetic order on increasing the oxygen 
xo parameter and, in particular, the antiferromagnetic ground state of the observed structure 
thus arise from the combined effect of an increased overlap of the magnetic orbitals and a sup- 
pression of the ferromagnetic local vanadium moment. Still, the importance of ferromagnetic 
correlations has been recently revealed by neutron scattering experiments, which suggested 
to characterize LiV204 as an itinerant d electron system close to weak ferromagnctism [ p2| . 
Similar conclusions have been drawn from the anomalous temperature dependence of the^Ji 
NMR relaxation rate |0, 12 . Finally, the very small values of the magnetic energy gain 



as compared to other transition metal oxides and the almost linear variation of this energy 
gain with the magnetic moment for the ferromagnetic states result from the minor role of a 
direct d-d exchange interaction as well as the small superexchange coupling due to 7r-type 
p-d overlap. Thus, without the geometric frustration, we would find an antiferromagnetic 
ground state, in which the well localized vanadium magnetic moments are coupled by weak 
superexchange interaction via the O 2p states. 

Nevertheless, the geometric frustration inherent in the spinel structure with only nearest 
neighbour interactions does prohibit long-range antiferromagnetic order as it leads to an 
energetical degeneracy of different antiferromagnetic configurations. Perhaps, in the states 
with large band width from the direct V-V overlap, the frustration of the magnetic interactions 
allows large amplitude spin fluctuations to occur. This could explain the T'^lnT behaviour 
of the specific heat observed at low temperatures, as well as the strong field dependence 
of the magnetic susceptibility below 30 K |g^. These spin fluctuations can result in large 
quasi-particle mass enhancements and thus mimic a heavy fermion state. Alternatively, in the 
absence of magnetic ordering, the states, in which the small overlap of the d with the O 2p 
orbitals occurs, could place the orbitals at the brink of a transition to a localization and the 
local spin fluctuations could be responsible for the observed heavy fermion behaviour. 

To summarize, the electronic properties of LiV204 as resulting from electronic structure 
calculations are dominated by a combination of broad and narrow bands near the Fermi energy, 
which originate from strong metal-metal and rather weak metal-oxygen overlap, respectively. 
While the direct V-V bonding does not contribute to magnetic ordering, the antiferromagnetic 
as well as the slightly unstable ferromagnetic state result from weak superexchange interaction. 
However, longe range magnetic order is suppressed by the geometric frustration inherent in 
the spinel structure; hence only large spin fluctuations survive, which provide a clue to an 
understanding of the observed specific heat and susceptibility data. The picture emerging 
from the present calculations is different from that usually used for / electron systems, where 
the hybridization between localized and the partially occupied itinerant states causes local spin 
fluctuations responsible for the HF behaviour, and when magnetic magnetic ordering occurs 
it is usually due to the RKKY interaction. 
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